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First study of B — > 7r semileptonic decay form factors using NRQCD 

H. Matsufuru a *, S. Hashimoto bt , K-I. Ishikawa a , T. Onogi a , and N. Yamada a 
a Department of Physics, Hiroshima University, Higashi-Hiroshima 739, Japan 
b High Energy Accelerator Research Organization(KEK), Tsukuba 305, Japan 

We present a quenched calculation of the form factors of the semileptonic weak decay B — » nW with 0(l/mq) 
NRQCD heavy quark and Wilson light quark on a 16 3 x 32 lattice at (3 = 5.8. The form factors are evaluated at 
six heavy quark masses, in the range of tuq ~ 1.5 — 8 GeV. I/tuq dependence of matrix elements are investigated 
and compared with HQET predictions. We observe clear signal for the form factors near q 2 ma x, even at the 
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1. Introduction 

Lattice study of B decay matrix elements is 
important for the determination of Cabibbo- 
Kobayashi-Maskawa matrix elements, and for in- 
vestigations of applicability of Heavy quark effec- 
tive theory ( HQET ) which is extensively applied 
to phenomenological studies. In this work, we cal- 
culate B — > 7T form factors using the heavy quark 
described by 0(l/m Q ) NRQCD and Wilson light 
quark jj]. 

2. Simulation 

Numerical simulation is performed on 120 con- 
figurations of a 16 3 x 32 lattice at (3 = 5.8 in the 
quenched approximation. Wilson light quark is 
employed with three values of k, 0.1570, 0.1585 
and 0.1600, which leads k c = 0.163461(69) and 
lattice cutoff a -1 = 1.714(63) GeV determined 
from p meson mass. Light quark field is normal- 
ized with the factor ^Jl — j|. 

The heavy quark is described with NRQCD in- 
cluding O(lmg) corrections which leads following 
evolution equation. 
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where A^ 2 ) denotes the lattice Laplacian and B 
is the chromomagnetic field. The stabilizing pa- 
rameter n should satisfy n > 3/2rriQ. 

For the heavy quark, eight values of mass 
and stabilizing parameter are used: (toq , n) = 
(5.0,1), (2.6,1), (2.1,1), (2.1,2), (1.5,2), (1.2,2), 
(1.2,3), and (0.9,2). m Q = 2.6 and 0.9 
roughly correspond to the b- and c-quark masses. 
The mean- field improvement || is applied to 
the heavy quark evolution equation with u = 
(!EW) 1/4 = 0.867994(13). 

The matrix elements are extracted from three 
point correlation functions, 



i(k—p) ■ x s 



C^(p,k;t f ,t s ,t l )=J2Yl e""e-* 

Xf X 3 

x<O fl (a: / )7 /1 t (aw)O ir (ti > 0)). 

We use 20 rotationally nonequivalent sets of (p, tt) 
with \p\, \k\ < a/3 • 27r/16. The source and the 
current operators are set on the time slices ti = 4 
and t s = 14 respectively. The matrix elements 
are extracted in the region tf = 23 — 28. 
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Numerical simulations were carried out on In- 
tel Paragon XP/S at INSAM (Institute for Nu- 
merical Simulation and Applied Mathematics) in 
Hiroshima University. 

We estimate the effect of perturbative correc- 
tions to the heavy quark self-energy and the cur- 
rent Q. In some cases, we use larger values of 
n for the perturbation than those in the Simula- 
tion: (mg.n) = (5.0,1), (2.6,2), (2.1,2), (1.5,3), 
(1.2,3), and (0.9,6). This is because of the 
singularities encountered in the perturbative ex- 
pressions for some set of (mg,n) with small n. 
The multiplicative part of the current renormal- 
ization constant is calculated with massless Wi- 
son quark for vanishing external momenta. Two 
scales q* — 7r / 'a and 1/a are considered to define 
the expansion parameter g v . 

3. Results 

It is useful to define following quantity. 
n(k)\V,\B(p)) 



W>*0 = 



y/2E„(k)y/2E B (jp) 



(3) 



This expression can be entirely composed of nu- 
merical results, without any assumption such as 
a dispersion relation. It is also convenient for a 
comparison with HQET predictions. According 
the heavy quark symmetry, for p 
constant value in the leading order of 1/toq: 
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where k\ — 2sin(fcj/2). Upper two of Figure 1 
show the results of V4 for p = k = and \% for 
|fc| = 27r/16 in the case of k = 0.1570. They are 
evaluated at three renormalization scales, mean- 
field tree, q* = n/a and 1/a. Both V4 and V% 
less depend on m B in comparison with fs case 
H . The spacial component of V is more affected 
by the perturbative corrections than the temporal 
one is. It is also predicted that 



V p (p = 0,k) 
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Figure 1. raj dependence of V for n = 0.1570. 
V4 (top figure), Vk (middle), and msVp (bottom) 
with three renormalization scales for each. 
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We extrapolate V p at finite p to p = linearly 
in p 2 to determine V p (p=0, k). V p (p = 0, fc = 0) 
multiplied by is also displayed in Figure 1. 
Contrary to the cases of V4 and Vk, 0(l/m_e) 
effect is significant for V p . 

The matrix elements are expressed in terms of 
two form factors, f° and / + : 



(7r(fc)|V r „|B(p)) = [p+k 



m ; 



-ml 



3 



m„=2.6, k=0.1570, Tree 
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Figure 2. Form factors for mq — 2.6 and K = 
0.1570 in the mean-field tree level. Solid lines 
correspond to the results of single pole fit. 



l -f(q 2 ). (7) 



To obtain the form factors from V^, it is neces- 
sary to assume some dispersion relations for n 
and B. For both of B and ir, we employ the 
form E(p) 2 = to 2 + J2i 4sin 2 (pi/2) to assure the 
on-shell conditions. Figure || shows the obtained 
form factors for mq — 2.6, k — 0.1570, in the 
mean-field tree level. It is observed that for larger 
mQ, / + increase rapidly near q max , which may 
be explained as the B* pole contribution. We fit 
the form factors to a single pole function, whose 
result is displayed in Figure || as the solid lines. 
The fit of / + gives the exchanged vector parti- 
cle mass mp i e — 2.979(49) which is comparable 
with tub* = 3.2502(62) determined from the two 
point correlator. The renormalization corrections 
change the shape of / + more than /°, because of 
asymmetric effect on temporal and spacial com- 
ponents. Even after the renormalization is incor- 
porated, the pole fit gives similar m po i e . 

We extrapolate the matrix elements linearly in 
X/k to the chiral limit. Alternative extrapolation 
using the form factors gives 7 % discrepancy of 
f°{q 2 max ) for m Q = 2.6 and 14 % for m Q = 0.9. 
At the zero recoil point, the prediction of heavy 
meson effective theory [|J advocates the linear ex- 
trapolation of the matrix element. The results at 
a£ q are much noisy for quantitative discussions. 
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Figure 3. Comparison of f°(q 2 max ) with fs/f-K in 
two scales, mean-field tree and q* = 1/a. 



Finally, we consider the implication of soft 
pion theorem 0,^). For the massless pion limit, 
f°(q 2 max) should equal to Jb/ fn- This relation is 
examined in Figure [|, using our result on fs de- 
termined with slightly different form of NRQCD 
p. Significant difference is observed in large ms 
region. Similar result is obtained in the work us- 
ing Fermilab action for the heavy quark J8| . The 
origin and physical meaning of this discrepancy 
remains as a future problem. 
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